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Abstract: In the midst of the rapid developments in electronic instruments and remote 
sensing technologies, airborne three-line array sensors and their applications are being 
widely promoted and plentiful research related to data processing and high precision 
geo-referencing technologies is under way. The exterior orientation parameters (EOPs), 
which are measured by the integrated positioning and orientation system (POS) of airborne 
three-line sensors, however, have inevitable systematic errors, so the level of precision of 
direct geo-referencing is not sufficiently accurate for surveying and mapping applications. 
Consequently, a few ground control points are necessary to refine the exterior orientation 
parameters, and this paper will discuss bundle block adjustment models based on the 
systematic error compensation and the orientation image, considering the principle of an 
image sensor and the characteristics of the integrated POS. Unlike the models available in 
the literature, which mainly use a quaternion to represent the rotation matrix of exterior 
orientation, three rotation angles are directly used in order to effectively model and 
eliminate the systematic errors of the POS observations. Very good experimental results 
have been achieved with several real datasets that verify the correctness and effectiveness 
of the proposed adjustment models. 

Keywords: airborne three-line array imagery; block adjustment; rotation angle; systematic 
error compensation; orientation image 
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1. Introduction 

The principle of three-line scanner imagery was first proposed by Hoffinan et al. [1]. The problem 
of the low vertical accuracy of adjustment, caused by the traditional small format images, was 
significantly improved by this method and has been successfully applied in the German MOMS02 
project with very good results [2-4]. In recent years, all of the high resolution Earth observation 
satellites, such as SPOTS [5,6], QuickBird [7,8], IKONOS [7-9], and more recently WorldView and 
GeoEye etc. , have utilized a linear push broom scanner to collect image data. Three-line scanners have 
been widely exploited not only in satellite photogrammetry, but also in aerial photogrammetry. In 
2000, LH Systems successfully appHed the three-line imagery technology in aerial photogrammetry 
and announced the airborne digital three-line scanner ADS40 [7,10,11]; and several airborne digital 
three-line scanners were released thereafter, such as Starlmager [12-14] and 3-DAS-l [14]. These 
airborne digital three-line scanners usually have a global positioning system (GPS) and inertial 
measurement unit (IMU) on board that can acquire the exterior orientation parameters (EOPs) with 
considerable precision and provide more observations for bundle adjustment. Under these conditions, 
the coordinates of ground objects can be precisely determined with only four control points distributed 
at the four comers, thereby providing a new solution for the automation of photogrammetry. The 
related processes such as image data preprocessing, triangulation [15-17] and fast generation of digital 
surface model (DSM) [6,18,19] have been thoroughly investigated. GPS/IMU aided processing of 
airborne three-line imagery becomes one of the central issues in aerial photogrammetry. 

Image data acquired by three-line scanners have a large base-height ratio, but the linear array push 
broom strategy makes the EOPs of scanning lines vary by time, and each scanning line thus has 
different EOPs. It is impossible to simultaneously solve the EOPs of all of the lines in bundle 
adjustment [1]. In order to reduce the unknowns, a proper sensor model must be established to 
represent the position and attitude of the sensor. Several sensor models, for example, the direct 
geo-referencing model [20-22], the piecewise polynomial model [1,23] and the orientation image 
model [1,13,24], have been proposed by researchers. The introduction of IMU has provided a great 
number of useful observations. Many researches to date have focused on using the cp-co-K rotation 
angle system defined by sequential rotations about the Y, X, and Z axes or a quaternion [15-17,24] to 
represent the attitude parameters. However, the attitude data output from IMU are usually based on the 
co-cp-K rotation angle system (or namely the roll, pitch and yaw system) defined by sequential rotations 
about the X, Y, and Z axes [25]. For better modeling the systematic error of IMU, this paper employs 
the co-cp-K rotation angle system and adopts both the systematic error compensation model and the 
orientation image model in the bundle adjustment. Four real datasets are used for our experiments, and 
the accuracies of both the systematic error compensation model and the orientation image model are 
discussed. The achieved results are also compared with the available results in the literature based on 
co-cp-K rotation angle system and quaternion. 
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2. Mathematical Model of Bundle Block Adjustment 

2.1. Basic Theory 

As previously mentioned, it is impossible to set the EOPs of all of the scanning lines as unknowns 
during bundle adjustment. Proper sensor models are necessary to approximate these unknown 
parameters. Generally, the direct geo-referencing model, systematic error compensation model, 
piecewise polynomial model, and orientation image model are often adopted. The standard output of 
IMU is based on the co-(p-K rotation angle system (or roll, pitch and yaw system) [24]. Therefore, only 
in the case where the recommended co-cp-K system is adopted in bundle adjustment can the identified 
systematic error in the IMU observations be eliminated via systematic error modeling. Detailed proof 
will be given in the next section. 

Regardless of which sensor model is chosen, the object point, the corresponding image point, and 
the photographic center should lie on the same line. Given the EOPs {Xsj, Ysj, Zsj, coj, (pj, kJ) of the j th 
scanning line, the instantaneous coUinearity equation of the push broom imaging system can be 
expressed as follows: 

fli {X -Xsj) + b^{Y- Ysj) + c^{Z- Z«,) 



^0 ~ f 



y-y^=-f 



{X - Xs) + (7 - Ys) + c,{Z- Zsj) 
[X - Xs,) + b^{Y- Ys) + c^{Z- Zs) 
a^{X- Xsj) + (7 - Ys) + c^{Z- Zs) 



(1) 



where (x, y) are the image plane coordinates of the image point; (xo, yo, f) are the interior elements; 
(Z, Y, Z) are coordinates of the object point; (X5,, Ysj, Zsj) are the translation parameters of the scanning 
line; and a„ Z>„ c, {i= 1, 2, 3) are the nine elements of the rotation matrix computed by the three angles 
{cOj, (pj, Kj) of the scanning line. 

2.2. Triangulation Based on the Systematic Error Compensation Model 

The positioning and orientation system (POS) integrated on the three-line digital scanner often has 
good precision and high stability. The post-processing modules will compensate for the offset related 
to the projection center and the misalignment related to the IMU main axes, but many experiments 
have shown that some residual systematic errors will remain in the corrected POS data [23]. It can 
be assumed to mainly consist of GPS and IMU systematic drift errors that vary by time. Taking the 
GPS antenna offset errors and the IMU primary axes misalignment errors into consideration, strip 
related constant and linear terms are used in this paper to represent the offset and drift errors of the 
GPS/IMU observations. The mathematical model to calculate the photographic center {Xs, Ys, Zs) with 
GPS observations {Xgps, Ygps, Zgps) is described as follows [26]: 
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where (m, v, w) are the remaining errors of the GPS antenna offset components; R is the rotation matrix 
of the angular elements {co, cp, k); {ax, ay, az) and {bx, by, bz) are the systematic offset and drift errors 
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of GPS observations of each strip; t is the imaging time of a certain scanning Hne; and to is the 
reference time that usually set to be the center time of a strip. 

Calculation of the angular elements involves the rotation matrix R/mi/ derived from IMU 
observations {coi, cpi, k/), the misalignment matrix Rmis derived from the remaining misalignment errors 
(com, (Pm, Km) of IMU primary axes, and the rotation matrix R derived from angular elements (co, (p, k) 
of the EOPs. These three rotation matrices are all based on the co-(p-K rotation angle system and have 
the following relationship [24]: 

R = Rm/sR/mc/ (3) 

Considering the strip offset parameters (a^u, a^, a^), and first order drift parameters (ba,, b^, b^) of the 
IMU observations, the systematic error compensated IMU observations (coj, cpi, Ki,) used in this paper is 
modeled as follows [26]: 
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where t is the imaging time of a certain scanning line; and is the reference time of a strip. 

The unknowns of the bundle adjustment include the remaining errors of the GPS antenna offset, 
the remaining misalignment errors of IMU primary axes, the strip offset, and the drift errors of the 
GPS and IMU observations, and the space coordinates of the object points. Given an area with m sfrips 
and n object points, the amount of unknowns is3 + 3 + 12xm+3 x n. After adjustment with the least 
squares principle, the exact values of the systematic errors can be obtained to refine the EOPs of all of 
the scanning lines. 

The proof of using the co-cp-K rotation angle system to model the systematic errors of IMU 
observations will be given by a numerical example. Suppose the IMU observations (<o/, (pi, ki), 
systematic errors (am, acp, Gk), and (ba, b^, b^) are coi = 0.050000, coi = 0.050000, kj = 1.550000; 

= l.Oe , = l.Oe , = l.Oe ; and = 5.0e , b^ = 5.0e , bK = 5.0e . According to the 
characteristics of IMU, the flying time of each strip is usually no longer than 20 min, so the maximum 
time offset of a scanner line is t - to = 600.00 seconds since to is usually set to be the center line time 
of a strip. 

According to the above supposed values, the systematic error compensated three angles (co,, cp,, k,,) 
can be calculated as follows: 
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(5) 



R/M(7 can be derived from the IMU observations (coi, (pi, kj), and liiMu can be computed by the 
systematic error compensated three angles (coi, cp,, k,,). 
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^0.0207688 -0.9985343 -0.0499792^ 
0.9984824 0.0232662 -0.0499167 
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0.0548713 -0.0529810 
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-0.0539738^ 
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0.9970868^ 



(6) 



Two sets of three angles {(pr, cot, Kt), and (qjir, cot, Kt) under the definition of the cp-co-K rotation angle 
system can be decomposed fi-om the above computed rotation matrix R/mi/ and R/m/, respectively: 



^r= 0.050062, cot= 0.049937, Kt= 1.547499 
(pr = 0.054079, c«'r = 0.05392 1 , /c, = 1 .55 1 082 



(7) 
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If the computed systematic errors are directly used as those under the definition of the (p-co-K 
rotation angle system, we can easily deduce the following three angles {(pc a>c, Kc): 



(8) 



Thus we can get the differences among the two sets of three angles: cp t ~ (p c ~ 0.000017, 
(Ot - (o'c = -0.000016, Kt - Kc = "0.000417. The focus to pixel size ratio of ADS40/80 sensors is 
always about 62.7/0.0065 = 9,630, so the maximum influence of pitch and roll errors caused by using 
different rotation angles is 0.000017 x 9,630 = 0.16 pixels, which usually can be neglected. However, 
the maximum influence of the yaw error at the image border is 6,000 x 0.000417 = 2.50 pixels 
considering the CCD length of ADS 40/80 sensor is 12,000 pixels, which is obviously unable to be 
neglected. This comparison shows that using the (p-co-K rotation angle system will cause errors of about 
2.50 pixels at the image border, because the strip-related constant and drift errors in IMU observations 
cannot be correctly treated under this angular system. Up to now, it is clear that modeling the 
systematic errors of IMU observations based on the co-cp-K rotation angle system is advantageous since 
the standard output of IMU device is under this angular system. 

2.3. Triangulation Based on the Orientation Image Model 

The basic principle of the orientation image model is to solve the three translation and three rotation 
elements of the EOPs at each orientation image by bundle adjustment. The EOPs of any scanning line 
between the two EOPs of the orientation images can be interpolated. First or third order Lagrange 
interpolation [24] is commonly used to interpolate the EOPs for airborne three-line imagery. The POS 
observations of an airborne three-line scanner have small random errors within a short time. Therefore, 
these observations can be used in the interpolation of EOPs. 

Given a certain ground point Pi , the corresponding image point {xij, yy) and the EOPs (Xsj, Ysj, Zsj, 
coj, (pj, Kj) of the corresponding scanning line j, the coUinearity relationship among them can be 
described with the EOPs of two adjacent orientation images k and A: + 1 as follows [24]: 
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Xy = ¥y {X. ,Yi,Zi, Xsk, Ysk, Zsk, co^,^^,/c^, Xsk + ,, Ysk + ,, Zsk + 1, ca^^^ ,(p^^^, k^^^ ) 

yy = G .. {X. ,Y^,Z^, Xsk, Ysk, Zs,, o)^,(p^,K^, Xsk . i, Ysk .,,Zs,.^ co^^^ ,(p^^^, k^^^ ) (9) 

where Fy and Gy represent the coUinearity equation of image coordinates {xy, yp), (Xf, Yi, Z,) are the 
spatial coordinates of point Pi, (Xsk, Ysk, Zsk, m, 9k, Kk) and (Xsk+i, Ysk+u Zsk+u m+h (Pk+i, Kk+i) are the 
EOPs of two adjacent orientation images k and k+ I. 

During the adjustment process, (Xsj, Ysj, Zsj, coj, (pj, kj) can be calculated by the EOPs of two 
adjacent orientation images and the correction values which can be derived by GPS/IMU observations. 
A variable Q is introduced to represent the relationship among the EOPs at time tj and the EOPs of the 
adjacent orientation images at time tk and tk+i- 

(10) 



'~t -t 



'■k+X '■k 

Thus {Xsj, Ysj, Zsj, coj, cpj, Kj) can be calculated by the following equations [24]: 

Xg. = CjXs, + (1 - )Xs, . , - AXsj 

Y, =C Ys, + (1 - C . )75a , 1 - AYs, 
Zg. = CjZsk + (1 - Cj )Zsk , 1 - AZsj 

cOj = Cjco, + (1 - Cj)co,,, - Am. 

^j = Cjrc^^{^-C;)K^,,-AK. 

where the correction vector AZ,, AYj, AZj, Acoj, Acpj, Akj, of the EOPs can be calculated by the EOPs at 
time j, k,k + I directly obtained from the GPS/IMU observations, as shown in Figure 1 [24]: 
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(12) 



The POS observations can be easily integrated into bundle adjustment with the orientation image 
model. The mathematical model is quite similar to the model described in the above section except that 
only the offset and drift parameters of the GPS/IMU observations should be taken into consideration. 
The unknowns of bundle adjustment now include the EOPs of the orientation images, the remaining 
GPS antenna bias components, the remaining IMU misalignment components, the offset and drift 
parameters of the GPS and IMU observations, and the coordinates of the ground objects. 

The existing strategies for bundle adjustment of three-line imagery typically adopt a quaternion [15,24] 
or Lagrange interpolation [13,18] to represent the attitude parameters. There is a constraint condition 
among the four elements of a quaternion. Usually, there are two solutions to deal with the constraint. 
The first solution is keeping the four unknowns and one additional constraint during bundle adjustment. 
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and the second solution is computing the fourth element with the other three elements calculated by 
bundle adjustment. 



Figure 1. Example of one orientation parameter over time [24]. 
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In this paper, bundle adjustment is realized based on the co-cp-K rotation angle system. The 
derivation of unknown coefficients related to translation elements is the same as that of the 
conventional method, whereas the derivation of unknown coefficients of IMU misalignment 
components {com, (Pm, km) and attitude parameters {cok, (pk, Kk), {cok+i, (pk+\, Kk+\) should be obtained by 
linearizing the coUinearity equations according to the m-cp-K rotation angle system. 

Given an block area with m strips, n unknown ground points and 1 orientation images, the amount 
of unknowns of bundle adjustment is3 + 3 + 12xm + 3xn + 6xl. The number of unknowns is 
extremely large during bundle adjustment, and thus the compression of normal equation and a fast 
matrix inversion algorithm are needed. EOPs of each scanning line can be interpolated by the corrected 
EOPs of orientation images. 

2. 4. Resolving the Interpolation Problem of Rotation Angles 

The range of angles in the three-line sensor data is usually between (-tt, tt), therefore, the periodic 
issues of angle interpolation must be solved when rotation angle is used to represent the attitude of the 
sensor. The values of roll co and pitch cp are definitely between {-nIA, nIA) in all cases since the optical 
axis of the image sensor is always pointing to the ground, but the yaw angle k is closely related to the 
fiight direction of the plane. For instance, when the plane flies in the direction of south to north or 
north to south, k is usually around ±n. The unconformity of the sign of yaw angles usually leads to 
wrong interpolated values. Actually, it is impossible for the plane to make a sharp turn that causes 
a big variation of k up to 180 degrees during a small period of time (the time interval between 
two orientation images). Hence, a simple but effective strategy by making the two angles have the 
same sign is proposed to resolve the interpolation problem. Given kf" , kf+i as the yaw angles of 
two orientation images at time k, and k+ \, the yaw angle A:f'^at time tj can be correctly interpolated by 
preprocessing of k'"" according to the following criteria: 
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jIMU _ rlMU 
tIMU _ jIMU 



+ 271, 
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7tl2 and kZ <-7r/2 



(13) 



The proposed mathematical model is simpler than that uses quaternion since the constraint 
condition among the four elements of a quaternion can be removed, which makes dealing with large 
dataset possible, although the computational complexity is increased. Moreover, this strategy is 
advantageous for the improvement of current bundle adjustment software that adopts rotation angles. 

3. Experiments and Analysis 

A bundle adjustment software program of airborne three-line scanner imagery named 
iBundle-AeroTLS® based on the proposed models was developed. In order to verify the correctness 
and effectiveness of the two models, ADS40 images collected from the Taigu and Pingyao test areas in 
China and the Waldkirch test area in Switzerland were used for experiments. The focal length and 
physical pixel size of ADS40 camera are 62.7 mm and 0.0065 mm, respectively. The adjustment 
results of the three test areas have been reported in papers [16,17] and are used for comparison in this 
paper. According to the characteristics of the ADS40 sensor [24], the time interval of the orientation 
image model was set at eight seconds during bundle adjustment. Additionally, the literature indicates 
that only four control points at four corners of the block are needed to obtain satisfactory results. Thus, 
the four control point strategy was adopted for all three experiments. 

The ADS40 images of the Taigu test area were captured in May 2007. The entire test area consists 
of 13 strips and covers an area of about 30 km . The ground sampling distance (GSD) is about 0.06 m 
while the flying height is about 600 m. A total of 105 ground control points (GCPs) were measured by 
differential GPS. Four control points on the four comers were used in the bundle adjustment and the 
other 101 points were used as check points. A total of 96,330 pass points were automatically matched 
by the photogrammetric software DPGrid developed by Wuhan University [27]. The bundle 
adjustment statistics for the two models are shown in Table 1, where RMSE is the root mean squared 
error of control and check point residuals, Mean and Max are the average and maximum of the 
residuals. The standard errors of unit weight in image space are 0.0042 mm and 0.0035 mm, respectively. 

Table 1. Error statistics of bundle adjustment of the Taigu test field (Unit: m). 



Point 
number 


Coordinates 


Systematic error 
compensation model 


Orientation image 
model 


Systematic error 
compensation and 
orientation image model 






RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


GCPs 
4 


X 


0.057 


0.019 


-0.073 


0.045 


0.013 


-0.061 


0.042 


0.013 


-0.058 


Y 


0.032 


0.002 


0.054 


0.031 


0.005 


0.055 


0.031 


0.005 


0.056 




Z 


0.050 


0.010 


-0.104 


0.044 


-0.012 


-0.095 


0.045 


0.011 


-0.091 


Check 


X 


0.045 


0.021 


0.103 


0.043 


0.015 


0.085 


0.041 


0.013 


0.085 


points 


Y 


0.039 


0.008 


-0.115 


0.035 


0.006 


-0.077 


0.035 


0.006 


-0.072 


101 


Z 


0.058 


-0.015 


-0.158 


0.046 


0.012 


0.115 


0.043 


0.012 


0.104 



The ADS40 images of the Pingyao test area were captured in October 2006. The entire test area 
consists of nine strips and covers an area of about 240 km . The GSD is about 0.2 m while the flying 
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height is about 2000 m. A total of 84 GCPs were measured by differential GPS. Four control points 
were also used in the bundle adjustment, and the other 80 points are used as check points. A total 
of 54,792 pass points were automatically matched by DPGrid. The bundle adjustment statistics for the 
two models are shown in Table 2. The standard errors of unit weight in image space are 0.0035 mm 
and 0.0030 mm, respectively. 



Table 2. Error statistics of bundle adjustment of Pingyao test field (Unit: m). 



Point 
number 


Coordinates 


Systematic error 
compensation model 


Orientation image 
model 


Systematic error 
compensation and 
orientation image model 






RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


GCPs 
4 


X 


0.105 


0.018 


0.206 


0.093 


0.012 


0.233 


0.091 


0.011 


0.206 


Y 


0.135 


-0.010 


0.286 


0.114 


0.020 


0.227 


0.110 


0.016 


0.201 




Z 


0.164 


0.040 


-0.323 


0.135 


0.025 


0.272 


0.123 


0.021 


0.237 


Check 


X 


0.129 


0.036 


0.328 


0.105 


-0.015 


-0.275 


0.100 


-0.015 


-0.275 


points 


Y 


0.148 


-0.025 


-0.332 


0.108 


0.018 


0.282 


0.098 


0.018 


0.252 


80 


Z 


0.185 


0.048 


0.487 


0.148 


0.022 


0.276 


0.143 


0.021 


0.271 



The ADS40 images of the Waldkirch test area with four east-west strips and two north-south strips 
were captured in May 2002. The ground coverage of this test area is about 80 km". The GSD of the 
imagery is 0.2 m while the flying height is about 2000 m and the number of measured GCPs is 30, 
respectively. All of the automatically matched 25,140 conjugate points among the six strips were used 
in the adjustment. Four control points at the four comers were used to perform the bundle adjustment, 
and the other 26 points were used as check points. The bundle adjustment statistics for the two models 
are shown in Table 3. The standard errors of unit weight in image space are 0.0038 mm and 
0.0035 mm, respectively. 



Table 3. Error statistics of bundle adjustment of Waldkirch test field (Unit: m). 



Point 
nmnber 


Coordinates 


Systematic error 
compensation model 


Orientation image model 


Systematic error 
compensation and 
orientation image model 






RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


RMSE 


Mean 


Max 


GCPs 


X 


0.124 


0.036 


0.241 


0.101 


-0.012 


0.215 


0.096 


-0.012 


0.203 


Y 


0.113 


0.051 


-0.213 


0.107 


-0.031 


-0.204 


0.112 


-0.022 


-0.182 


4 


Z 


0.142 


0.044 


0.223 


0.125 


0.021 


0.215 


0.121 


0.014 


0.202 


Check 


X 


0.087 


0.048 


0.184 


0.087 


0.037 


0.185 


0.085 


0.032 


0.159 


points 


Y 


0.088 


0.045 


0.216 


0.078 


0.018 


0.217 


0.076 


0.015 


0.210 


26 


Z 


0.162 


-0.070 


-0.279 


0.149 


-0.021 


-0.275 


0.145 


-0.015 


-0.253 



In the above three experiments, the differences of the estimated misalignments of GPS and IMU 
obtained by the two approaches are well below 0.02 m. However, the position and attitude elements 
of EOPs of the orientation image model have changed a few centimeters and several arc seconds 
respectively when compared with those of the systematic error compensation model. As shown in 
Tables 1-3, the RMSE values of the check points in planimetry were smaller than 0.7 GSD, and the 
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maximum residuals ranged from 0.7 GSD to 2.0 GSD after bundle adjustment with the systematic 
error compensation model in the three test areas. The RMSE values in the height direction were 
slightly better than 1.0 GSD, and the maximum residuals were about 1.3~2.5 GSD. However, some 
residual systematic errors remained in both the planimetry and height. After bundle adjustment with 
the orientation image model, the systematic components of the check point residuals were significantly 
reduced, for example the absolute mean value of height residuals of check points in Table 3 decreases 
fi-om 0.070 m to 0.021 m. The RMSEs in the height direction substantially decreased to smaller than 
0.75 GSD in the three test areas. Moreover, the maximum residuals in both the horizontal and vertical 
directions were smaller than 2.0 GSD, and the remained residuals show random distribution. 

The results of combining both the systematic error compensation model and orientation image 
model were also listed in the three tables. It was shown that about 5 to 10 percent improvements were 
achieved when compared with using the orientation image model only. 

The available results of bundle adjustment for the three test areas in the literature [16,17] were 
listed in Table 4. By comparison among the results listed in the four tables, we can see that the 
accuracies of bundle adjustment with the systematic error compensation model for the three test areas 
were equivalent to those with the commercial software ORIMA. However, the accuracies of bundle 
adjustment with the orientation image model were better than those calculated by ORIMA. For 
example, the RMSEs of height residuals were 0.046 m and 0.148 m in Tables 1-2 respectively, while 
they were 0.06 m and 0.19 m in Table 4, which means that the results of our approach with the 
orientation image model are up to 20 percent higher than those of the literature. These satisfying 
results verify the correctness and feasibility of the proposed two adjustment models based on co-cp-K 
rotation angle system. 



Table 4. Adjustment results from published reference papers [16,17] (Unit: m). 



Test field 


Point number 


RMSE 


Maximum residual 


X Y Z 


X Y Z 


Taigu 

PINGYAO 

Waldkirch 


4 GCPs 101 Checkpoints 
4 GCPs 84 Check points 
4 GCPs 26 Check points 


0.07 0.05 0.06 
0.15 0.19 0.19 
0.174 0.174 0.304 


0.12 0.10 0.13 
0.45 0.43 0.45 
N/A N/A N/A 



To fully demonstrate the advantages of using the <x)-(p-K system against the quaternions, a large 
project with 41 ADS40 strips that covered more than 7,000 km was used for experiments. The project 
was located in Henan Province (China). The 41 strips' images with 0.45 m GSD were acquired by 
three flight missions. The amount of automatically matched conjugate points by DPGrid was more 
than 300,000 for this test dataset. There were 41 field measured ground points, in which 13 points were 
used as GCPs and the other 28 points as check points in the bundle adjustment as shown in Figure 2. 
The size of absolute orientation tiles (*.odf), which contained the exterior orientation elements of each 
scanner line, was about 874 MB. There were totally 4,312 orientation images by using the orientation 
image model with 8 s intervals. The band width of unknowns was 2,196 after optimization. 
Consequently, the needed memory to store the coefficient matrix of the normal equation was about 
434 MB (4,312x6x2,196x8/1,024/1,024 = 434 MB). Our proposed approach and developed software 
was capable to deal with the large dataset as a whole project with both systematic error compensation 
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model and orientation image model and both gave stable and similar bundle adjustment results. The 
result of bundle adjustment with orientation image model was listed in Table 5. The RMSE of height 
residuals was about 1.0 GSD, and the maximum residual in height was about 2.0 GSD. There was no 
stereo view hardware available while performing this experiment, so it was not possible to precisely 
measure the image points under stereo view environment. Although the height accuracy was 
comparatively lower than those of the other three test fields because the image coordinates of GCPs 
and check points were not precisely measured, it still qualified for the requirements of topographic 
mapping at scale 1:5,000 in China. 



Table 5. Error statistics of bundle adjustment of Henan province test field (Unit: m). 



Number of points 


RMSE 


Mean residual 


Maximum residual 


X Y Z 


X Y Z 


X Y Z 


13 GCPs 
28 Check points 


0.450 0.520 0.464 
0.296 0.476 0.408 


0.043 0.034 0.065 
0.012 0.024 0.106 


-1.126 1.052 -0.857 
0.734 1.006 0.924 



Figure 2. Distribution of fiight lines, GCPs and check points of an ADS40 project with 
4 1 strips. ■ means control point, o means check point. 



5 



f 

4. Conclusions 

For airborne three-line scanner imagery, aerial triangulation is a prerequisite of photogrammetric 
product generation. Rotation angles were adopted in this paper to express the attitudes of the sensor, 
and mathematical models of systematic error compensation and orientation image were proposed for 
aerial triangulation. Using the rotation angles against quaternions has the advantages of eliminating the 
systematic errors effectively and dealing with large datasets. 



o 
o 
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The results of three test blocks demonstrate that the RMSE values of the check points are smaller 
than 1.0 GSD in both the planimetry and the height after bundle adjustment with the systematic error 
compensation model. However, there are still some systematic residuals for the check points since the 
systematic errors of the GPS/IMU observations could not be fully compensated. The accuracies of the 
check points in both the planimetry and the height were substantially improved by bundle adjustment 
with the orientation image model; the RMSE values in the height direction were refined in particular, 
fi-om 1.0 GSD to 0.75 GSD. The systematic error of the check point residuals was also compensated, 
which demonstrates the reasonability and superiority of the orientation image model. Generally, 
5 to 10 percent improvements were achieved by the combination of the two models when compared 
with using the orientation image model only. Moreover, dense conjugate points with very good 
quality automatically matched by DPGrid software also contribute to improve the accuracy of 
bundle adjustment. 
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